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Abstract. I investigate the correlation between spatial and flow anisotropy in determining the ellip- 
\ tic flow and azimuthal dependence of the HBT correlation radii in non-central nuclear collisions. It 

£NJ ■ is shown that the correlation radii are in most cases dominantly sensitive to the anisotropy in space. 

' In case of V2, the correlation depends strongly on particle species. A procedure for disentangling the 

spatial and the flow anisotropy is proposed. 
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In non-central nuclear collisions at RHIC energies, the resulting fireball can exhibit 
anisotropy in both spatial shape and transverse expansion velocity profile. They both 
influence the measured "elliptic flow" coefficient V2 III]. A question arises: how are 
they correlated in the determination of V2, i.e., which combinations of spatial and flow 
anisotropy lead to the same elliptic flow? 

On the other hand, dependence of HBT correlation radii on the azimuthal angle is also 
shaped by the two mentioned anisotropics. Therefore, the same question can be asked: 
how do the spatial anisotropy and transverse expansion flow anisotropy combine in the 
(j) -dependence of correlation radii? 

An analogical situation appears in determining the slopes of single-particle pj- 
spectra. It is well known that they are determined by temperature and transverse ex- 
pansion velocity and that it is impossible to disentangle these two quantities from a 
single measured spectrum. There is, however, also the -dependence of HBT radii in 
which the correlation of temperature and transverse flow is qualitatively different from 
that in the determination of spectra. Temperature and transverse flow velocity then can 
be unambiguously measured from analysing both spectra and HBT radii. 

A similar solution shall be sought here: can we disentangle spatial and flow anisotropy 
in non-central collisions by analysing both V2 and the azimuthally sensitive HBT radii? 

Note that several statements have been made in literature which are related to this 
programme. In [2] the STAR collaboration concluded that it was impossible to determine 
spatial anisotropy just from the measurement of V2 and a conjecture was made that HBT 
analysis would be able to gain such result. Two qualitatively different final states resulted 
from hydrodynamic simulations by Heinz and Kolb [3] and the authors demonstrated 
the possibility to distinguish these states by HBT interferometry. Here I report on a 
systematic study of the interplay between spatial and flow anisotropy in framework of 
generalisations of the blast-wave model. 



FIGURE 1. The two different models for transverse expansion velocity used here. 



AN AZIMUTHALLY ANISOTROPIC BLAST- WAVE MODEL 

Instead of fully describing the used model I will just focus on those features which are 
important for this work and refer the reader to literature for more detailed discussion 
10,111]. Suffice it to say that the fireball is thermalised with a temperature T and exhibits 
longitudinally boost-invariant expansion. Its transverse profile is ellipsoidal and the 
emission function is 




oc 0(1 -r), (1) 



where R x and R y are the two transverse radii, in and out-of-plane, respectively. They can 
be parametrised with the help of a spatial anisotropy parameter a 

R x = aR, R y = -. (2) 
a 

Thus an out-of-plane elongated source is characterised by a < 1, whereas for an in-plane 
elongated source we have a > 1 . 

The transverse expansion velocity also depends on the azimuthal angle. The velocity 
is given as 

v ± =tanhp(f,0). (3) 

We shall have a closer look at two models which differ in the azimuthal variation of 
the velocity. In Model 1 [5] the velocity is always perpendicular to a surface given by 
f = const. This direction together with the reaction plane defines the azimuthal angle (j)/,, 
as illustrated in Figure[T] The transverse rapidity 

p(f,0)=fpo(l+P2Cos(2fo)), (4) 



where the parameter po measures the radial flow and p2 is the flow anisotropy parameter. 
As the velocity is perpendicular to the surface of the fireball, this model resembles the 



expansion profile early in the fireball evolution: the direction of velocity coincides with 
acceleration which in turn is given by the pressure gradient. 

In Model 2 the transverse expansion velocity is directed radially and varies with the 
usual azimuthal angle, which is denoted as S here 

p(f,0) =fp o (l+p 2 cos(20 s )). (5) 



THE ELLIPTIC FLOW 

Recall that v 2 is defined as the second Fourier coefficient of the azimuthal dependence 
of spectrum 
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It can be calculated in the two used models and the result reads ||4|] 

y _ fodrrf^d<l> cos(20)7(0)gi( a )/ 2 (fr) 
V2 fid77fS*d+J(4)K 1 (a)I (b) ' 

where the arguments of the Bessel functions are a = mj coshp(f, 0)/T and b = 
pr sinhp(r, (j))/T. The only difference between the two models appears in the Jacobian 

m 

Model 1: J((j)) = (a 2 cos 2 (j> +a~ 2 sm 2 (j>) , (8a) 
Model 2: 7(0) = (a~ 2 cos 2 + a 2 sin 2 0) . (8b) 

From these relations it is obvious that the two models lead to the same v 2 if they are 
related by transformation a — > a ~ 1 . In other words, one in-plane and another out-of-plane 
source give the same v 2 . This is an analytic illustration of the claim that it is impossible 
to determine even the qualitative type of spatial anisotropy just from measurement ofv 2 . 

Now we can look on the correlation between flow and spatial anisotropy and study it 
only for Model 1, since results for the other Model are obtained simply by substitution 
a —* aT x . In Figure we see that the correlation between a and p 2 strongly depends 
on the particle species. Hence, here is a strategy for determining both a and p2~. first 
determine the temperature and radial flow coefficient po from azimuthally integrated 
spectra. Their dependence on azimuthal anisotropics was shown to be small o]. Then 
measure v 2 for at least two particle species and obtain a and p 2 . Of course, this procedure 
assumes that we know which model to use for the analysis. This leaves an open question 
which is to be answered by correlation measurement. 



AZIMUTHALLY SENSITIVE HBT 



In non-central collisions, the HBT correlation radii can be measured as a function of 
the azimuthal angle 0. We shall focus mainly on the two transverse radii R a and R s and 
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FIGURE 2. Elliptic flow V2 calculated with Model 1 for pions (upper row) and protons (lower row) 
with pt = 0.2GeV/c (left column) and 0.5 GeV/c (right column). The used values T = 100 MeV and 
Po = 0.88. Thickest contour lines show where v'2 vanishes; consecutive lines correspond to steps by 0.02. 

decompose their azimuthal angle dependence as 

R 2 o(<i>) = < o + 2< 2 cos20 + ... (9a) 
^ 2 (4>) = ^ 2 o + 2^ 2 ;2 cos20 + .... (9b) 

The individual terms of these decompositions are obtained as various combinations of 
space-time variances taken with the emission function |@, 0]. Because we are rather 
interested in the oscillation of the radii and not so much in their absolute size, we shall 
look at the normalised oscillation amplitudes R^/Rfo El 1 - They are sensitive to a and 
p2, but less sensitive to R and po- 

From Figure|3]we conclude that the azimuthal oscillations of the HBT correlation radii 
are mainly shaped by the spatial anisotropy parameter a. Dependence on flow anisotropy 
is weaker, with the only exception of at high K t in Model 2 which is determined 
mainly by flow. This confirms the statement that the azimuthal dependence of correlation 
radii follows mainly the spatial anisotropy, especially at low K t . This has been shown 
here in framework of two models. It would be natural to expect this behaviour to be 
valid in general. It can be spoilt by very strong flow gradients which differ by much in 
in-plane and out-of-plane directions. A questions arises, however, whether large enough 
difference of the flow gradients is realistic. 



1 In fact, Retiere and Lisa realised in 1 5] that because it also includes time contributions, R 2 is not a good 
normalisation quantity and so used /?•? to normalise all of R 2 2 , R 2 s2 ^ an d R^i 2- This is not done here. 
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FIGURE 3. Normalised oscillation amplitudes ^ 2 /^«o (upper row) and R^/^sO ower rows) as 
functions of a and p 2 calculated for T = 0.1 GeV, p = 0.88, R = 9.41 fm, To = 9fm/c, At = 1 fm/c 
(see |4J] for definitions of all parameters). Left columns show results for Kj = 0.3GeV/c, right columns 
correspond to Kj = 0.9 GeV / c. Thickest lines show where the second-order oscillation terms vanish, other 
contours are set in steps of 0. 1 . 



The two investigated models exhibit similar dependence on the spatial anisotropy 
parameter a when focusing on the oscillation of HBT radii. Recall, however, that they 
were related by transformation a — > aT x when reproducing the same vj- Therefore, two 
different models which both reproduce v% measurement will behave differently when 
fitting the azimuthal dependence of HBT radii. This is illustrated in Figure El Both 
models used in this figure fit measured the V2 for pions and protons well. However, 
while Model 1 reproduces the RHIC data qualitatively well, Model 2 leads to the phase 
of oscillation just opposite to data [8]. 

Thus we conclude that among the two models used in this study, Model 1 seems to 
correspond to RHIC data, whereas Model 2 is clearly ruled out. This does not disqualify 
it, however, from future applications at the LHC where possibly longer lived fireballs 
could be produced which will develop a different transverse flow pattern. 



CONCLUSIONS 

It has been demonstrated analytically that one cannot disentangle spatial and flow 
anisotropy of the fireball just from a measurement of vi. I also demonstrated that, at 
least for two classes of models, the azimuthal dependence of correlation radii reflects 
the type of spatial anisotropy the source actually exhibits. 

Thus I can propose the following (schematic) procedure for disentangling a and 
p2'. first measure the azimuthal dependence of HBT radii and determine the spatial 
anisotropy a. Then, with that a try to reproduce V2 for more species. Since for different 
species a and p2 are correlated in different ways, this should lead to unique pair of the 
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FIGURE 4. Comparison of the dependence of and R% on azimuthal angle (j) with the data measured 
by STAR collaboration (Au+Au at 200 AGeV, centrality class 20-30%) The curves and data points 
correspond from top to bottom to Kj = 0.2, 0.3, 0.4, 0.52GeV/c. Parameters of the models are: T = 
0.12GeV, p = 0.99, p 2 = 0.035, R = 9.41 fm, T = 5.02 fm/c, At = 2.9fm/c, and a = 0.946 (Model 1) 
or a = 1.057 (Model 2). Both models were tuned to fit v 2 {pt) ffl- 



anisotropy parameters. 
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